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I propose an alternative Hamiltonian field theoretical model for a light quark 
condensate that is compatible with QCD in the deep infrared. Key elec- 
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^ ■ I. INTRODUCTION 



^ . The aim of this paper is to reconsider again the old pre-QCD Nambu-Jona-Lasinio(NJL) 

ideas, but this time in the light of the new insights gained in the last 10 years. We shall make 
appropriate use of contemporary methods and ideas on boson and fermion condensates that 
proved to be so useful in both nuclear physics and condensed matter physics. The phenomena 
associated with broken chirai-flavour symmetries and chiral anomalies are indeed to be 
ultimately related to the physical nature and structure of a presumed quark-gluon quantum 
vacuum. 

We assume the existence of a u,d fermion condensate and work from there, using math- 
ematical methods more or less common to all fermionic systems at zero temperature. Our 
basic degrees of freedom are chiral QCDu,d,s-quarks but only in colour singlet combinations. 
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Their couplings to the leptons and gammas are assumed to be primarily as prescribed by 
the Standard Model (SM). 

Updating the NJL idea, we design an effective Hamiltonian to be worked out under 
the assumption that it has a stable Dirac-Hartree-Fock-Bogoliubov (DHFB) state as its 
(approximate) ground state . Necessary renormalizations are carried out using experimental 
meson masses and electroweak data. Applications and detailed numerical results will be 
published elsewhere. 

II. THE RELEVANT EFFECTIVE HAMILTONIAN 

We begin by making the standard symmetry assumptions [1,2, 3]. The internal symmetry 
that phenomenology suggests as the most relevant to the analysis of the lightest quark sector 
is the chiral-ffavour symmetry of the form 

G : SU NL (g) SU NR (g)U 1L (g)U 1R (1) 

with N=2,3. The chiral-flavour symmetry G undergoes various quantum mechanical sym- 
metry breakings [2,3] that reduces it down to flavour groups 

G =>- H : SU NL+R ® U 1L+R (2) 

Thus G is a symmetry of the Hamiltonian (or equations of motion) whereas H is a 
symmetry of the quantum vacuum. There are two well-understood kinds of symmetry 
breakings that achieve this, both traceable to the quantum vacuum [2, 3, 4]. The first kind is 
due to the existence of a specific classically definable Landau-like vacuum long range order 
(LRO)[5], to be defined more precisely in the following sections. It breaks the chiral SUn 
(SUn) part of the algebra (1), leading to the emergence of the right set of Goldstones(viz. 3 
and 8). The other kind of quantum mechanical symmetry breaking removes the axial U\l-r 
symmetry of the Hamiltonian, as a consequence of the axial anomaly. The latter can be 
interpreted as simply a consequence of the physical existence of the Dirac vacuum[4]. This 
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implies that there are no chiral doublets in the hadronic phase of QCD in which we live. 
This information is the input used in designing the effective Hamiltonian of this paper. 

We choose our basic degrees of freedom to be included in this effective Hamiltonian: 
these are (we use Bjorken and Drell notations, definitions and conventions [6]) 

(u L , d L ) © (u R , d R ) © (s L , s R ) (3) 

and their antiparticles.We assume input mass matrices: 



^ m u 



Mi 



On 



m d 
m, 

Thus the simplest non-trivial relevant effective Hamiltonian can be assumed to be 



V 



(4) 



H 



St 



U + H + V« + V< 2 > (5) 



Hp 



E / d 3 xq nR (x)(-tfV)q nR (x)) + L^R (6) 



V (1) = E [d 3 xq nR (x)M 0n q nR (x))+L^R (7) 

... „J 



n=cu,cd,cs 



V {2) = Yl J d3 ^QnL(x)q nR (x)q n/R (x)q n/L (x) +L R (8) 

X n,n'=cu,cd,cs 



9VA 



E / d 3 xq nL (x)YqnL(x)qn>L(xhuq n iL(x) + L R (9) 



8ttA2 , 

X n,n'=cu,cd,cs 

The q n L,R are chiral field operators [6]. The necessary inputs for this model are: 

(i) a real flavour diagonal 3x3 "input mass matrix" M 0n ; 

(ii) Two real independent effective dimensionless couplings gsp and gvA] 

(iii) two fundamental mass scales, provided by an explicit A x ~ lGeV (related to the 
scale where massless quarks become massive quarks) and an implicit Aqcd ~ 300MeV (a 
parameter basically fixing the overall size of physical hadrons). In the context of a non-chiral 
quark model this region would give the medium range qq potential. 



III. THE RELEVANT LRO 



Let us define the Landau-like long range order (LRO) parameter appropriate to an "u,d- 
quark condensate" by assuming the existence of a robust spinless, colourless and flavourless 
non- vanishing scalar LRO [2,3] 
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E < 0\(q nL (x)q nR (x))\0 > +L ^ R (10) 



8ttA2 

X n=cu,cd,cs 

Translation invariance ensures independence on space coordinates. This defines g^p.The no- 
tation |0 > as used here should be explained. In a many-body context the symbol |0>would 
usually mean that one is referring to a certain state vector in Hilbert space representing the 
true ground state of the many-body system. In the context of this paper, such a true ground 
state of QCD-quark-gluon coupled fields is of course not only unknown but it is also irrel- 
evant. This notation is nevertheless adopted here, but |0 > merely defines a "no particle 
state" which is just a tautology for "normal operator products" [7]. It has therefore nothing 
to do with any true physical ground state. We shall sometimes refer to it rather loosely as 
the " quark vaccum". 

We work exclusively with approximate Heisenberg operators that play the role of "physi- 
cal states" [7]. It is established that only colour singlets would qualify as such. We try to find 
stationary solutions to the Heisenberg equations of motion for these "physical states". Thus 
our "no particle state" is nothing but a tautology for normal ordering DHFB quasiparticle 
operators. 

In order to include this assumption in the effective Hamiltonian we begin by making a 
straightforward Boguliubov transformation using the Nambu-Gorkov representation for the 
u,d QCD-quarks[5]: 

\ / 
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(11) 



sin ip nh x(p) cos (fnhx(p) 
- cos (f nhX (p) sin <p nhX (p) 
where A = L,R are chiralities and h are helicities. The "Bogoliubov angles (p n hx(p)" serve 
as adjustable variational parameters. By separating out the bilinear from the non-bilinear 
terms we find that 



U 8t (u,d,s) = U^ + l^ + V' (12) 



where 



H o = E Jd 3 p:(q nR (p) q nL (p) 



( 



\p\ -A 

[~^n -\P\ ) 



QnL (P) 



(13) 



and V represents the remaining (quadrilinear) terms. The Bogoliubov angles are chosen so 
that Hq is fully diagonalized: 

K = -K+ E E + (14) 

n=cu,cd,cs h=±7£ 

sin 2 ipip) = \{l- \p\/E u4 (p))) n = u,d (15) 
E u4 (p) = v^+Kd + Af (16) 



E s (p) = J\p\* + ml (17) 



where by definition 



b nh (p)\0 >= = d n - h (-p)\0 > (18) 



Note the independence of the Bogoliubov angles on helicities / chiralities that follow from the 
definition of the LRO. We shall refer to this as the Dirac-Hartree-Fock-Bogoliubov(DHFB) 
static approximation. 

IV. RENORMALIZATIONS 



We shall have to renormalize the above theory, selecting experimental data on pions and 
the etas for doing so[l]. These mesons, just like any physical mesons, are here considered to 
be just complex poles of the physical S-matrix amplitudes. 

The electromagnetic sector of this model is represented by the Hamiltonian 



H ew — H / 7 + V ew (19) 

and will be treated perturbatively.H / 7 stands for free leptons and gammas. The Y ew is 
given by the SM[2,3]. 

(i) Consider the 7r ± main decay channel in the rest frame , defining the decay constant 

/tt± 

A( tt ± — > // ± + i/^)) =< 0| J wc «fc|7r ± (M w ±) > (20) 

where 

\^{M^ ± ) >= |tt+ >= ^EE / ^ ± (p) *6^^(-p)|0 > (21) 



< 0\J weak \i^(M n ±) > | = \f n± \JJ%±- (22) 



2(2tt) 



we find the condition 



M v ± \f w ±\ r d 3 p 



I 72^S*^=(P)*«» 2 V(P)I (23) 



V 3(2vr) 3 4cos^ 'i (2tt) 
(ii) Consider the 7r ± charge radius, defined through 

< vr + (p 2 )|JL(0)k + (Pi) >= G w ((p! -P2) 2 )^— ====== (24) 

where 

G w ((pi - P2) 2 ) = 1 + \ < rl > (px - p 2 ) 2 (25) 

A simple estimate of the theoretical charge radius can be obtained by making a reasonable 
ansatz for the (common) normalized internal wavefunction of n ± : 

*«±(p;x,y) = -J=N(x,y)eM-\(p- I^QCD? (26) 



where x,y are dimensionless variational parameters. So 

.j >= A 2^(wyAoCT) (27) 

*F 2 (x,y; A X /A QCD ) 1 J 
6 



/oo 
duexp(-(A x /A QCD ) 2 u 2 )(x + y U y (28) 
-x/y 

(iii) Next, consider the main decay channel of the 7r°at rest: 

A ^ _^ 7 (fc) +7 (-fc) =< 7 (£), 7 (-£)|V em |7r° > (29) 

k°>=iE E E/^ o (p1^(p)^(-p)|o> (so) 

V i<! c q=cu,cd h J 

But assuming the initial pseudoscalar at rest, we find that the absorptive part of the (dom- 
inant) underlying Lorentz and gauge invariant qq annihilation amplitude a (with massive 
quarks) is : 

k 2 Im a(7r°(^^ q + q) — ► l{k) + -f(-k)) ~ 

~ (A/A;) 2 (1 - (A/kfyh l n[ i _ (i _ (A/A;) 2 ) */(l + (1 - (A/A;) 2 )^] (31) 

.As A, A; — ► it can be shown that the lhs of (31) tends to 5(k 2 ) [8], contrary to naive 
expectations. The deep reason for this is the existence of the fixed anomaly pole at k — 0, 
a necessary feature if the Ui em gauge invariance is to be maintained. This pole lives below 
the physical threshold( at about k = 2 A), which can be related to the physical mass of the 
7T° through the Gell-Mann-Oakes-Renner formula [2,3]. The existence of this pole could be 
guessed from (31) as the rise of k 2 lma as k descends from infinity towards the physical 
threshold. This provides another condition on our parameters. 

(iv) In order to calculate the value of the 770 — i]' anomaly [2, 3] we shall have to include 
further interaction terms from (12). We shall consider the effect only to leading order and 
ignore all isospin mixings. Let us define 

fa >= -4EEE / d 3 p^ v MKf ih (p) d ifd-p))\° > 1 = x > 2 ( 32 ) 

V 6 c h fi 

with fi = u,d and f2 = s. 

So by diagonalizing the 2x2 matrix 
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Mij =< r}i\H st (u, d, s) - U \r}j > (33) 

we find that the eigenvectors (in the approximation of keeping only the qq pair exchange 
diagrams) are 

\r] >= cos#|??i > +sin6>|77 2 > (34) 

\r]' >= sm9\r]i > — cos9\r] 2 > (35) 
where the mixing angle is given by 

tantf = (M v o - M n )/M 12 = M 21 /(M % - M 22 ) (36) 

M u = 2E u4 +\F 1 \ 2 M 22 = 2E S + \F 2 \ 2 M 12 = F*F 2 + F 1 F* = M 21 (37) 

I r°° I r°° 

F^-^J^ d 3 p* m (p)M^P) F ^ = ^J d 3 p^ m (p)A 2 (M s ,p) (38) 

A 1 (A,p) = l2A/J\iP\+A* (39) 

A 2 (M s ,p) = M 8 /y/\p*\ + M? (40) 
The angle is thus part of our renormalized parameters. 

V. CONCLUSIONS 

We presented an Hamiltonian field theory in a DHFB-static approximation as a comple- 
ment/alternative to the conventional theory of pions and etas. However, our theory (based 
on an updated version of the NJL field theory[9]) can be easily extended and has a much 
broader scope. It can also easily be linked to non-chiral quark models and QCD in the 
deep infrared. All renormalizations procedures are build in the theory itself and have of 
course only meaning in the context of this theory, as every renormalization scheme does in 
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its own context. The issue of confinement, though not essential to our case, is bypassed in 
a natural way. Detailed numerical fits to experimenta data in order to get the renormalized 
parameters will be published elsewhere. 
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